Oscillatory optical response of amorphous plasmonic nanoparticle arrays 
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The optical response of metallic nanoparticle arrays is dominated by localized surface plasmon 
excitations and is the sum of individual particle contributions modified by inter-particle coupling 
depending on specific array geometry. Here we scrutinize how experimentally measured properties 
of large scale (30 mm 2 ) amorphous Au nanodisk arrays stem from single particle properties and their 
interaction. They give rise to a distinct oscillatory behavior of the plasmon peak position, full- width 
at half-maximum, and extinction efficiency which depends on the minimum particle center-to-center 
(CC) distance. 



PACS numbers: 78.67.-n, 41.20.-q, 42.25.Dd, 78.40.Pg 
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Strong coupling of light to metal nanoparticles via lo- 
calized surface plasmons is one reason for the wide ex- 
ploitation of nanosized metallic entities^ - — For many tar- 
geted uses of nanoplasmonic systems a key question is 
whether to operate with individual metallic structures^ 
or to use ensembles in the form of periodio^ or random 
arrays on a support 4 The optical properties of nanoplas- 
monic arrays, both periodic and fully random, stem from 
the optical response of individual particles. The array 
modifies these single-particle spectra, sometimes quite 
considerably, via inter-particle coupling that depends on 
the exact array geometry. Thus, array design generally 
is an additional handle for tuning plasmonic response, 
together with particle size, geometry, and materials. 

Here we scrutinize experimentally and theoretically a 
novel oscillatory behavior of the optical response from a 
specific type of nanoparticle array, somewhere between 
perfectly periodic and fully random, that we refer to as 
an amorphous array. This particular type of nanopar- 
ticle arrangement on a surface exhibits short range dis- 
tance order, while, at long distances, it is completely ran- 
dom. Furthermore, it can be quite easily fabricated on 
large areas (wafer scale), using bottom- up self-assembly 
based nanofabrication techniques like hole-mask colloidal 
lithography, 9 making it a first choice for many large- 
scale devices and applications^^ - — Our finding shows, 
in contrast to the generally accepted opinion, that amor- 
phous arrays exhibit distinct properties of interacting 
particles even if their density is low. 

For our experiments, we fabricated large area ar- 
rays of gold nanodisks with engineered randomness us- 
ing an electron beam-lithography (EBL) nanofabrication 
schemed Circular areas of roughly 30 mm 2 were pat- 
terned for each considered center-to-center (CC) distance 
with minimum imposed CC ranging from 2.5 to 7 in units 
of particle diameter D (C). The size distribution for the 
disks is very narrow as can be seen in Fig. [T]and the his- 
tograms in Ref£ and, consequently, basically eliminates 
inhomogeneous broadening. We fabricated large arrays 
(30 mm 2 ) with a very large number of particles to elimi- 
nate the effect that for the same set of global parameters 
(minimum CC, D, thickness 20 nm, illumination condi- 
tions fixed) small samples could correspond to slightly 
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FIG. 1: (color online) Experimentally measured extinction 
spectra of gold nanodisks with an engineered randomness. 
Resonance position, peak value, and linewidth all show a non- 
monotonic dependence on lattice parameter C being minimum 
center-to-center distance in units of disk diameter D: (a) D = 
160 nm, (b) 260 nm. (c) and (d) show SEM images of the 
amorphous arrays for C — 3 and 4, respectively. Notice the 
random distribution of perfectly defined particles. 



different array realizations. 

In the first two panels of Fig. [1] we present experi- 
mentally measured extinction efficiency spectra near the 
resonance to illustrate the sensitivity of the peak position 
and extinction efficiency per particle at peak to the CC 
value C for D = 160 nm in (a) and D = 260 nm in (b). 
In a first rough analysis for D = 160 nm we see, that for 
C = 2.5 extinction is maximal at 838 nm, then it redshifts 
for C = 3 and 4 to 855 and 867 nm, respectively, and then 
undergoes a blueshift to 837 nm (C = 5) and 815 nm for 
C = 6 - suggesting an oscillatory behavior of the peak 
position. A similar trend is also seen when tracing the 
peak amplitude (extinction efficiency) and the linewidth. 

An efficient way to model arrays of plasmonic nanopar- 
ticles is by a coupled dipole approximation in which each 
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disk is modeled by an induced point dipole coupled to 
an external electromagnetic fields In this framework, 
the particle properties are described by a polarizability a 
determined by the material, geometry, and surrounding 
mediumi 14 ' 15 In the quasistatic regime the polarizabil- 
ity a qs is proportional to V(e m - e s )/(e s + L(e m ~ e s )), 
where e m and e s are the permittivities of the metal parti- 
cle and surrounding medium, respectively, V is the par- 
ticle volume, and L is a shape depolarization factor. 
Dynamic depolarization and radiative damping are ac- 
counted for by introducing the modified long wavelength 
approximatio n 15 ' 16 1/a = l/a qs — |ifc 3 — where k is 
the wave number of exciting light of wavelength A and a 
is a length associated with the particle geometry.— 

For an infinite periodic array, where the particles are 
interacting, the system of coupled equations is solved 
by assuming that the polarization of each particle is 
the same and thus a becomes an effective polarizabil- 
ity that takes into account inter-particle interactions via 
a retarded dipole sumi£ ' 17 i 18 However, with a gradual in- 
crease of disorder, the narrow peak characteristic for a 
periodic array disappears and the response turns into an 
inhomogeneously broadened plasmon resonance^ 

For an amorphous array we can, as an ensemble av- 
erage, define an effective polarization a*. One way of 
analyzing the inter-particle contributions to a* is to av- 
erage over many realizations of amorphous arrays (i.e. 
dipole sums). Here, however, we describe a model in 
which the average particle is surrounded by a continu- 
ous film of dipoles with surface densities determined by 
the pair correlation function Q(r, C), where r is the radial 
distance from the considered particle. One can think of 
this as an average of an infinite number of different real- 
izations of amorphous arrays centered around a specified 
particle placed at a particular specified point. In a sense 
this approach is reminiscent of the coherent-potential ap- 
proximation for a random distribution of particles on a 
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square arrays 

For the average particle in an amorphous array we 
carry out the same procedure of solving the discrete 
dipole equations as ir>ii where the retarded dipole sum 
(discrete particles) is replaced by a retarded dipole inte- 
gral (continuous film with hole) 
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where the exponential term multiplied by the expression 
in the square brackets (e ikr [■■■]) describes the retarded 
dipole-dipole interaction, g(r,C) — crQ(r,C) is the pair 
correlation function Q(r, C) multiplied by the particle sur- 
face density tr = uq£~ 2 , £ cc = CD, and o~o is a surface 
packing parameter. The integration is over the whole 2D 
(r, 6) space with the exception of an inner circle smaller 



than £ cc . Performing the angular average yields the av- 
erage, effective polarizabilityii£ 
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where 
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The function e tkr (k 2 + (l — ikr)/r 2 ) consists of two parts: 
the first (e %kr k 2 ) comes from the far-field dipole radiation 
and its value oscillates, while the second (e lkr (1—ikr) / r 2 ) 
corresponds to intermediate and near-fields and its value 
has a well-defined limit for r — ¥ 00. However, these obser- 
vations are only strictly valid for a well behaved function 
G(r,C) at infinity. 

To perform the integration for S in Eq. [3] we require 
an expression for G(p, C), where p = r/D is a normalized 
radius. We obtain Q(p,C) by finding a function which 
fits well (R 2 ~ 1) to pair correlation data calculated 
from random distributions of particles generated with 
the random sequential adsorption algorithm, which was 
also used to calculate particle positions for the fabricated 
arrays*^ The selected function consists of two parts - a 
constant and a varying one 
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It is chosen because its product with functions describing 
dipole fields is relatively easy to compute and its R 2 = 
0.99 (fitting parameters given in22). Slight differences 
between this function and the one shown by Hinrichsen et 
al£^ occur only for particles at close distances. However, 
as we show later, a qualitative description is insensitive 
to the exact expression describing the short range order. 

Equation 01 while easily integratable when multiplied 
by the expression for dipole radiation does not lend itself 
to an easy exposition of the main physics taking place. 
Therefore a qualitative analysis is carried out first, before 
performing the full calculation. We do this by keeping 
the hard-core part (unity) and omitting the second, os- 
cillating term (multiplied by the sine function). Thus, the 
simplified Q reduces to a Heaviside step function 0(p — C) 
that describes a fully random array with a removed circle 
of radius C around the average particle. The simplifica- 
tion still describes satisfyingly the most important prop- 
erty of the array, namely the short-range order defined 
by the minimal allowed CC distance for the analyzed 
particle, and does not alter the main physical processes 
occurring within the array. 

The function Q = Q(p — C) is derived from Eq. H]by 
setting the pair correlation function to unity. To calcu- 
late the far-field term, which oscillates around a mean 



value, we modify it by adding to the exponent the term 
— er, which makes the expression nak 2 e lkr ~ Er dr 
well-defined. This represents a case when the array is 
illuminated by a very broad Gaussian beam, i.e. disks 
very far away from the center of the beam feel a dimin- 
ished intensity of the electric field, but the decay is slow 
enough that the average treatment of the disks holds. 

The first integral in lim e ^o equals ■nokie %M ^ and the 
second irae lMcc /£ cc . Thus, for the simplified case of Q = 
Q(p — C), the retarded dipole integral becomes 
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We substitute this result into the effective polarizabil- 
ity a* (Eq. [2J and rewrite the right hand side contain- 
ing the substituted expressions into a Lorentzian form to 
easily identify the peak position and full-width at half- 
maximum (FWHM). We employ, as an illustrative test 
case, a generic metal sphere with a Drude dielectric func- 
tion e(uj) = 1 — ujp/ (u)(uj + ij)). Using the modified long 
wavelength approximatio n 1 ^ 16 we get 
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where we have introduced dimensionless variables uj = 

2 

— , 7 = ■ SL , s = ^ with = -£ the Mie resonance fre- 
quency, a coupling strength q = it<7q(R/£ cc ) 3 (q is max- 
imum f for ao — 1), two functions f(k£ cc ) and g(k£ cc ): 
f(x) = cos x — x sin a;, g(x) — sin x + x cos x, and c is the 
speed of light. 

From the denominator of Eq. ^ we can read off 
the resonance frequency and the FWHM. In the non- 
interacting case (q = 0) we have a mode at 1 
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with a linewidth ^^ 



For an amor- 
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phous array consisting of such particles the resonance 
is modified by inter-particle coupling via the retarded 
dipole integral and Eq. ((SJ shows that we have a reso- 
nance at 



fli/n N = \/i - qf(k N £ cc ) 



(7) 



in terms of the non-interacting J7n and a linewidth of 

fi/f N = 1 + 2(1 + s 2 )g N5 (fc N 4c), (8) 

where Qn is the individual particle quality factor. We 
see from Eqs. © and flSJ that due to the quality fac- 
tor we expect the randomness to show up much more 
in the FWHM than in the resonance frequency. This is 
clearly seen in Fig. [5] which shows peak position (red) 
and FWHM (black) based on Eqs. © and ©. The 
peak position follows a sinusoidal line with a period de- 
termined by the ratio of the minimum particle-particle 
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FIG. 2: (color online) The variation of the optical charac- 
teristics of the particle arrays in Fig. [TJ can be understood 
from a simple generic Drude model. The resonance frequency 
(red) and FWHM (black) of the array using a hard core pair- 
correlation function with minimum center-center distance £ cc 
is shown. We measure the center-center distance £ cc in units 
of the resonance wavelength An- The resonance frequency as 
well as the linewidth are in units of the bare particle proper- 
ties. Notice how the interference between the particles causes 
both resonance frequency and FWHM to oscillate and the 
oscillations of the FWHM are larger than for the resonance 
frequency. With decreased coupling the curves settle to the 
particle values. 



distance to the resonance wavelength An of a single par- 
ticle. Notice, that the oscillations are governed by the low 
cut-off £ cc . The real part of the dipole interaction term 
shifts the resonance position towards higher frequency 
when 9?(<S) = / < 0, while for 5i(<S) > it induces a 
red shift. The largest shifts occur when the array is rel- 
atively dense, cf. large coupling strength q. For large 
inter-particle spacing the interference between the disks 
vanishes since q is proportional to £~ 3 . 

The single particle resonance linewidth T\ is modified 
by the imaginary part of S = g which introduces a 
modulation of the linewidth equal to (1 + s 2 )Q^g(kfq£ cc ). 
Similar to the peak position, the linewidth is a decay- 
ing oscillatory function that tends to the single particle 
resonance width for infinitely diluted amorphous arrays. 
When g < for even-numbered half-periods the reso- 
nance linewidth is smaller than the single particle one, 
however, it does not go to zeroi 15 ' 23 " — Note, that in this 
analysis we have kept the minimum center-center dis- 
tance larger than of the order of two diameters so that 
higher order multipoles which are not present in the the- 
ory should be negligible. 

In the above qualitative picture, in which we dropped 
the oscillatory term in the pair correlation function, we 
have shown that the oscillations of the optical cross sec- 
tions of amorphous arrays are the result of interference 
between the incident field driving a particle and the scat- 
tered fields originating from the other particles in the ar- 
ray. We use now the full expression for Q (Eq. to 
analyze experimental data, shown in Fig. [TJ obtained 
from extinction measurements on nanofabricated amor- 
phous arrays of nanodisks. To model the properties of a 
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experiment model 
■ D= 480 nm — 
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FIG. 3: (color online) Experimental results compared to our 
predictions based on the full pair-correlation function for the 
extinction at peak value, the peak position and linewidth as 
a function of the minimum particle-particle distance in units 
of particle diameter. The characteristics of the spectra (ex- 
tinction value, position, and linewidth) oscillate as a result 
of radiative coupling between the plasmonic particles within 
the array and are a function of the minimum CC distance, 
(a) shows peak extinction values for a disk in the amorphous 
array for the different disk diameters, which are normalized 
to the peak extinction values of their respective single par- 
ticles; (b)-(d) show extinction peak position (red left y-axis) 
and linewidth (blue right y-axis) for diameters (b) 160, (c) 
260, and (d) 480 nm as a function of the CC distance. The 
horizontal thin dash-dotted lines indicate asymptotes of the 
peak position and linewidth of the amorphous arrays. 



single disk in the array, we adjust an oblate spheroidal 
polarizability so that the single disk (for a hypotheti- 
cal case of an infinitely diluted array) resonance position 
and linewidth correspond to the asymptotes of the ex- 
perimental arrays for very large CC. 

First, we address the extinction cross sections per par- 
ticle in the arrays. Figure [3ji presents the measured and 
calculated extinction per disk in the array, C c , normal- 
ized to the extinction of a single disk. Clearly, extinc- 
tion exhibits strong oscillatory CC distance dependence 
for all three measured particle sizes. The minimum-to- 
maximum difference is about 40% for small CC distances 
and agrees very well with the model calculations. 

The measured experimental peak positions (diamonds) 
and full-width at half-maximum (FWHM, circles) values 
are shown in Fig. [3j>d. The solid red and dashed blue 
lines representing peak position and FWHM values, re- 
spectively, are calculated according to the scheme out- 
lined above. The agreement between the theoretical and 
experimental data is very good, in particular in view of 
our relatively simple theoretical treatment of the particle 
interaction in the array. The period and phase shift of the 
measured and calculated oscillations of the peak position 
and FWHM are consistent and show a pronounced CC 
distance dependence. Notably, as the most extreme case, 



the experimentally measured FWHM for the D = 260 
nm disk varies between ca. 0.3 and 0.5 eV at a rather 
moderate change of the peak position (1.15 to 1.25 eV). 

We now briefly address trends in the amplitude of the 
oscillations that physically originate from radiative cou- 
pling as a function of nanodisk size in the amorphous 
array. It is known from theory and experiment that the 
scattering efficiency, especially compared with absorp- 
tion (i.e. the scattering/absorption branching ratio), in- 
creases with particle size£& Thus, it is to be expected 
that the oscillations are most pronounced for the largest 
particles, where the radiative coupling is strong, and de- 
crease in amplitude as the nanodisk diameter decreases. 
This is precisely the trend that can be seen in Fig. [3] 
(for D = 480 nm the maximum peak position oscilla- 
tion amplitude (AE) normalized to peak position (Eq) 
is AE/E = 0.11, for D = 260 nm AE/E = 0.09, and 
for D = 160 nm AE/E = 0.06). Consequently, the de- 
crease of the oscillatory amplitude is expected to continue 
as the nanoparticle diameter further decreases. This can 
be clearly seen in Fig. [2 where a 60 nm Drude sphere 
was considered. 

The observed discrepancies between the calculated val- 
ues for the peak position and, in particular, the FWHM 
are the result of several factors, among which are in- 
homogeneity of the fabricated disks (this effect is, how- 
ever, small), the fact that disks with large diameter-to- 
thickness ratios are not perfectly described by one dipole, 
and the estimation of the long-range interference term in 
the dipole integral. Furthermore, higher order terms also 
start to become important. Another issue to be noted 
here is the illumination of the array in the modeling, i.e. 
the assumption of a very broad Gaussian beam incident 
onto an infinite amorphous array (the array is in focus, 
so the phase of the incident beam is uniform). The latter 
assumption is then used to calculate the far-field term. 
However, for a finite array this may not be fully correct. 
Far-field radiation is proportional to e lkr and its definite 
integral (present in the retarded dipole integral) oscil- 
lates, so the value of the sum for a finite array depends 
on the exact relation between the array and the illumi- 
nation. Thus, our result should be viewed as an average 
value which lies between a minimum and maximum and 
may be a little off for some arrays. 

In summary, we have shown experimentally and ex- 
plained, using a dipolar model, an oscillatory optical re- 
sponse of amorphous plasmonic nanoparticle arrays that 
depends on the minimum allowed particle-particle sepa- 
ration. Optical spectra of amorphous arrays, while stem- 
ming from those of single particles, exhibit a strong in- 
fluence of intra-array radiative coupling of the plasmonic 
disks that results in oscillation of the extinction, its po- 
sition and linewidth. 
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